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Abstract 

We study the formulation of bond-orientational order in an ar- 
bitrary two dimensional geometry. We find that bond-orientational 
order is properly formulated within the framework of differential ge- 
ometry with torsion. The torsion reflects the intrinsic frustration for 
two-dimensional crystals with arbitrary geometry. Within a Debye- 
Huckel approximation, torsion may be identified as the density of dis- 
locations. Changes in the geometry of the system cause a reorgani- 
zation of the torsion density that preserves bond-orientational order. 
As a byproduct, we are able to derive several identities involving the 
topology, defect density and geometric invariants such as Gaussian 
curvature. The formalism is used to derive the general free energy for 
a 2D sample of arbitrary geometry, both in the crystalline and hexatic 
phases. Applications to conical and spherical geometries are briefly 
addressed. 
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1 Introduction 



According to the KTNHY theory @> 0> two-dimensional melting in the 
plane is a two stage defect-driven process involving continuous crystalline- 
to-hexatic and hexatic-to-fluid transitions. The intermediate hexatic phase 
is characterized by quasi-long-range bond orientational order and positional 
disorder. In a wide variety of settings one may expect to encounter geome- 
tries more general than the plane. In the theory of membranes the geometry 
itself fluctuates. External forces may act to bend the geometry to some fixed 
curved surface. This raises the challenging problem of generalizing the estab- 
lished KTNHY theory to substrates with some arbitrary geometry. In this 
paper we discuss a complete geometrical formulation of bond-orientational 
order [|J (hexatic and crystalline) that provides a framework in which the 
interaction of defects, geometry and topology are easily formulated. 

To begin, consider a hypothetical flat monolayer displaying bond-orien- 
tational order at zero temperature. Bonds may be represented as six vectors 
at each point of the plane, with a relative angle of |. The six vectors at any 
given point of the plane are parallel to those at any other point. In other 
words, if we know the vectors at some point we can reconstruct the vectors 
at any other point by parallelism. Imagine now adiabatically deforming the 
plane of the monolayer to an arbitrary curved surface. Intuitively we would 
expect the bond-orientational order to be stable to this deformation. The 
curvature of the monolayer, however, now implies an associated Gaussian 
curvature. As well known from differential geometry there is no intrinsic 
notion of parallelism on curved surfaces. To define parallel transport we 
require a rule specified via a connection. In general parallel transport of 
a vector between two points will depend on the path. For our problem 
this means that the bond angle at an arbitrary point on the surface has no 
path-independent meaning. For the standard choice of connection (the Levi- 
Civita connection) the bond angle at some point reached by two distinct 
paths from a reference point will depend on the total Gaussian curvature 
enclosed by the paths. This sensitive dependence on the path suggests that 
bond-orientational order is incompatible with any curved geometry. 

This mathematical argument, however, contradicts our intuition that 
bond-orientational order should be stable to deformation of the surface. 
This suggests that a more physical connection exists for defining bond- 
orientational order on curved surfaces. To see this note that bond-orien- 
tational order very generally implies the existence of a frame of six bonds 
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Figure 1: The result of parallel transport of V from P to P' a) phase with 
no disclinations: the result is path-independent b) phase with a few isolated 
disclinations, marked as grey spots: the result may differ by an angle of | 
but bond-orientational order is preserved c) isotropic fluid phase: parallel 
transport is completely ambiguous. 

forming gle at each point of the sample. Given this data we can define 

a path-independent parallel transport in the following way. Take a vector V 
at point P. This forms a certain angle with the frame. Now we can transport 
V to the point P 1 by demanding that V forms the same angle to the given 
frame at P' . By construction this physical parallel transport does not depend 
on the path and is therefore unambiguous []. 

Clearly the connection that corresponds to the parallel transport defined 
above must be different from the standard Levi-Civita connection. We will 
argue below that the mathematical realization of this physical connection 
appropriate for defining bond-orientational order involves a new degree of 
freedom, a geometrical quantity called torsion. The difference angle between 
two parallel transported vectors will depend on the enclosed curvature of the 
new connection with torsion. This connection may be tuned to be curvature- 
free so that parallel transport is path- independent, as desired. 

lr The fact that this parallel transport is unambiguous can be made more explicitly if, 
before starting the experiment, we mark all the bonds pointing, say, positively in the 
x-direction, which can be done unambiguously in the absence of free disclinations. 
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Although torsion may appear at this stage as merely a mathematical trick, 
a number of authors J7j have identified torsion as necessary for describing a 
system with a high density of dislocations, such as in the hexatic phase. 
This is certainly the case if the density of dislocations is well approximated 
by a continuous density, corresponding to the Debye-Huckel approximation. 
We will see more generally, however, that torsion accounts for the intrinsic 
frustration of a 2d configuration on a curved surface. 

This curvature of the connection with torsion also has a clear physical 
interpretation as being proportional to the density of disclinations. Raising 
the temperature slightly will excite disclinations. Although the zero curva- 
ture condition is violated, it is such that parallel transport is ambiguous only 
up to an angle 7r/3, and the hexatic order is preserved. As the temperature 
is raised further the growing density of free disclinations will eventually melt 
the hexatic to a fluid. The high density of disclinations in this phase may 
be represented by a continuous disclination density. In this case the paral- 
lel transported bond angle is completely ambiguous, which is equivalent to 
saying that bond angles are meaningless in the isotropic fluid phase. 

The organization of the paper is as follows: in the next section we intro- 
duce the necessary background in differential geometry, which will also serve 
to fix our notation. In section [3] we provide the physical interpretation of 
torsion and derive some general relations involving the geometry, the defect 
density and the torsion. In section ^ we reformulate the case of a simple 
flat monolayer in terms of the new formalism: this is easily generalized to 
an arbitrary geometry. We end the paper with some conclusions. In the ap- 
pendix we treat the example of a cone displaying hexatic order as an explicit 
application of our formalism. 

2 Differential geometry background 

In this section we briefly emphasize those concepts in differential geometry 
that play a major role in the subsequent analysis. For further details we refer 
the reader to the literature 0. 

2.1 Basic differential geometry 

A D-dimensional surface embedded in <i-dimensional Euclidean space is de- 
scribed by c/-dimensional functions i? a (x) , a = l,---,d of D coordinates 
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x M , [i = 1, • • • D, with D = 2 being the physically interesting case. 

The fact that the surface is embedded in H d provides a natural metric, 

ds 2 = g^dx^ <g> dx v , g^ u = d^R ■ d u R . (1) 

Introducing the vielbeins, e a fJi , it is possible to rewrite the previous equation 

as 

ds 2 = 5 ah $ a <g> $ b , (2) 
where $ a = e a fl dx^. It is easy to check that the tangent vectors 

e^e^A, a = i,..., jD , (3) 

where e [i a e a v = g^, define a basis of orthonormal vectors. The dot product 
of any two vectors u = u a e a and v = v a e a in tangent space is u • v = 5 a bU a v b , 
as if the metric were flat. This formalism is called the non-coordinate basis. 
Hereafter we use the first letters of the Latin alphabet when we work in a 
non-coordinate basis, and Greek indices for the coordinate basis. 

The space of all r-dimensional forms of a manifold M Vt r (M) plays an 
important role. The metric allows one to define a natural isomorphism, the 
Hodge-star operator 

Q r (M) A n D ~ r (M) , (4) 
defined as the linear map which in a non-coordinate acts as 

1 

(D 

where e 12 "' D = 1. We may define an inner product in Q r (M). If to — 
±LU ai ... ar ^ A • • • $ ar and v — ±is ai ... ar $ ai A ■ ■ ■ r d a,r 

(lu\v)= [ ujA*u= ( uj ai ... ar u ai "' ar \f\g\dx 1 ■ ■ ■ dx D , (6) 
Jm Jm v 

the integral being over the whole manifold M. 

The adjoint d^ of the exterior derivative is defined as the linear operator 
satisfying 

(da\P) = {a\S(5) . (7) 

It is an easy exercise to check that d^ = (_i) Z)r + D + 1 * d*. The Laplacian is 
the operator 

A = dS + Sd . (8) 
Let us recall that the Laplacian is independent of the connection. 



* (tf ai A tf a2 A • • • A $ ar ) = — -r-€ ai - ar ar+1 ... aD ('d ar + 1 A • • • A f D ), (5) 
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2.2 Connections 



A connection on a manifold M specifies how tensors are transported along a 
curve. It allows one to define a covariant derivative V M on a vector field U, 

v,u» = d,u» + r; a u° , (9) 

where r^ CT are the connection coefficients. The parallel transport of a vector 
along a curve &• is fixed by demanding that is be covariantly constant: 

V c y = ^P-VvV* = . (10) 
at 

It is convenient to define the connection coefficients in a non-coordinate basis: 

V a e b = T c ab e c , (11) 

and a connection one-form via 

cA = r a cb r . (12) 

We consider only metric connections that preserve the norm of a vector under 
parallel transport, a property called metric compatibility. This property 
constrains the connection one-form to satisfy 

oo\ = -u\ , (13) 

which in D = 2 further simplifies to cu a b = — e a ftf2, Q = £l a $ a - Connection 
coefficients do not transform as a tensor but it is possible to construct two 
tensorial objects out of them, the torsion and the curvature. For the sake of 
completeness, we write them acting on vector fields X, Y, 

T(X,Y) = V X Y-V Y X-[X,Y] 

R(X,Y) = VxVy - Vy Vx - V[x,V] (14) 

Expressions turn out simpler in a non-coordinate basis. Defining torsion and 
curvature 2-forms by 

T a = -T a b J b A$ c 

R\ = ^#A^, (15) 
5 



they are related to the connection one-form by Cartan's structure equations, 
which read for the case D = 2, 



T a = cftf° - 1M tf 6 e° 6 (16) 
-e a b R a b = dfl . 

Recall that the curvature 2-form is locally an exact form. The expressions in 
the coordinate basis follow easily from those in the non-coordinate basis from 
the tensorial nature of curvature and torsion, namely T^ up = e tJl a e b u e c p T a i )C 
and R^ 1 vpa 6^ a 6 u & p6 tjR bed- 



2.3 Some global considerations 

The global structure of the manifold places important constraints on the 
differential geometry of the surface. From Eq. [16] the curvature 2-form is 
locally an exact form, but that is not the case globally. For a manifold 
without boundaries, integrating over the whole manifold we have 

/ \t a b R a b = 2^Ind a = 2n X , (17) 
Jm I a 

where Ind are the indices of the vector field associated with Q and x is the 
Euler characteristic of the manifold. If R is the Riemannian curvature this 
is the usual Gauss-Bonnet theorem. 



3 The geometry of bond-orientational order 

We now consider a two-dimensional sample exhibiting bond-orientational or- 
der and forming an arbitrary shape in d-dimensional space, specified by the 
embedding R a (x). Within the sample, distances are measured according to 
the metric Eq. [L| 

The main feature of connections with non-vanishing torsion (non-sym- 
metric connections) is that infinitesimal reference vectors fail to close when 
parallel transported along each other, as illustrated in Fig. |[ If a crystal is 
forced to lie on an arbitrary geometry, the atoms in the crystal cannot be all 
equally separated by a distance a. This implies that parallel transport along 
the bonds joining the nearest-neighbors will fail to close, as shown in Fig. 0. 
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Thus we see that torsion is related to the intrinsic frustration experienced by 
a 2D crystal on an arbitrary surface. 

In the hexatic phase there is a similar interpretation of torsion, except 
that there will be a large number of dislocations. This is specified by a 
Burgers vector density 

b(x) = -^£bf5(x 4 -x) (18) 

where hf is the microscopic Burgers vectors with origin at point Xj. In 
a highly dislocated medium, the discreteness of the Burgers vector density 
b(x) may be approximated by a continuous vector density, the so-called 
Debye-Huckel approximation. 




Figure 2: Parallel Transport in the presence/absence of Torsion. 

The density h(x) is, in fact, directly related to the torsion invariant ap- 
pearing in Eq. [15], as noted by several authors 0. The torsion two- form 
is 

T a = v a ti l A $ 2 = v a tt M , (19) 

where Qm = ylol^jdx 1 A dx 2 is the volume form. Since v a is a vector it can 
be written in a coordinate basis = e^ a v a . Its geometric interpretation is 
illustrated in Fig. 0. Let us assume non-zero torsion at point A. Consider 
the vector AC and parallel transport it along vector AB, the resulting vector 
is BD'. If we parallel transport vector AB along AC we get CD" . Using the 
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parallel transport Eq. [II] and Eq. ^, we may write to lowest order 

WD 7 ^ = {T£ u - T^}AB' T AC U = v» X Area of the region (ABD'D"C). 

(20) 

If v = then D' = D" = D and the parallelogram closes. 

It is clear from the previous geometrical argument that torsion is related 
to the continuous density of Burgers vector by 

b(x) = v(x) , (21) 

or T a = b a QM, in a non-coordinate basis. A highly dislocated media phys- 
ically is represented geometrically by a surface endowed with a connection 
with torsion 0. 

Clearly then torsion is a necessary ingredient in a geometrical description 
of bond-orientational order, as it allows a connection form from which the 
bond angle is constructed by parallel transport. This interpretation of torsion 
is completely general. In the Debye-Huckel limit we have the additional 
identification of the torsion with the dislocation density. 

In the absence of free disclinations we pointed out that this result should 
be independent of the path chosen. This in turn translates into the mathe- 
matical requirement of vanishing curvature 2-form 

R a b = . (22) 

This is a local condition, not always possible to fulfill globally, as for example 
in the case of spherical topology. Anyway, Eq. [22] is fulfilled everywhere for 
the case of a sample having the topology of a plane. 

In this regime, the Cartan structure equations Eq. [T6| read 

T a = d& a - n A $ b e a b , (23) 
= dVL , 

The last equation implies Q = d9, where 8(x) is a 0-form, a function. Let 
us interpret its physical meaning. The equation of parallel transport Eq. |10| 
reads 

dya n rfc V V b - (24) 
where c^(t) parametrizes a curve in the surface joining the points P and P' . 

V ± (P') = ViiPy-^c^" , (25) 
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In the hexatic phase with no free disclinations, Eq. ^3] implies 

V ± (P') = V ± {P)e~^ p,) , (26) 



where V± = V 1 ± iV 2 and 0(P) = 0. This parallel transport id manifestly 
path independent and 9 clearly is the bond angle at point P' . Dragging P' 
along the entire surface and performing the parallel transport Eq. ^6| we can 
unambiguously construct the bond angle at any point from the knowledge of 
the bond angle at P. We have succeeded in giving a precise mathematical 
formulation of the heuristic zero temperature experiment performed in the 
introduction. 



3.1 Introduction of free disclinations 

The sample may also contain some free disclinations, either as a result of 
thermal fluctuations or because topological constraints force them to ap- 
pear. In either case, the parallel transport can only be defined modulo an 
angle of |. The system responds to the raising of the temperature by gener- 
ating curvature, but only via defects that preserve overall hexatic order, the 
disclinations. In this way, the curvature is directly related to the density of 
disclinations @. 

From these arguments, it is clear that a finite density of N free disclina- 
tions is represented as 

\e a b R\ = s(x.)n M 
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7T 1 N 



a(x) = x^E<^(x-x,), (27) 

6 y\9\ j=i 

where qi may be +1 (a five-fold disclination) or —1 (a seven-fold disclination). 
Higher integer charges are also possible although rarely need to be considered. 
If the surface is closed, we can integrate the previous equation over the whole 
manifold. With the aid of the Gauss-Bonnet theorem Eq. |I7] we find 

E q 3 = 6* (28) 

3=1 

This shows that, apart from the torus with x = 0, closed manifolds always 
contain a certain number of free disclinations. 
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With free disclinations, the Cartan structure equations (Eq. |T^) become 

T a = d& a - VI A $ b e a b , 

dtt = s(x)tt M . (29) 

Away from the actual location of the disclinations Q = dd, but now it is 
generally impossible for 6 to be defined as a continuous function. It is easy 
to check that a vectors parallel transported from point P to P' following two 
different paths c(t) and f(i) paths (see Fig. |XJ) differ by an angle ip 

^ = ^E«/> ( 30 ) 

where j runs over all charges within the region enclosed by the two paths. 
The ambiguity is a multiple of | and as expected hexatic order is preserved. 
Mathematically we have a connected manifold with a Z$ holonomy. Discli- 
nations may be regarded as orbifolds. 
Generally, we will write 

n = do + n sing , (31) 

explicitly displaying the bond angle 9 and a a singular part, obviously cor- 
responding to a vortex contribution. This is nothing but the procedure of 
separating a regular part and a singular part in the standard treatment of 
the XY model [|. 

There is more information encoded in the Cartan structure equations 
(Eq. |16|). First of all, the Levi-Civita connection is obtained imposing the 
vanishing of torsion, 

dti a — Vt L A $ b e a b = , (32) 

where Q L is the connection form of the Levi-Civita connection. We have the 
important relation 

dQ L = ^e\Rf = K(x)n M , (33) 

where R^ b is the standard Riemann tensor two-form and K(x) is the Gaus- 
sian curvature of the surface. Using this relation in Eq. [29| 

Q L A $ b e a b - Q A $ h e a b = T a . (34) 

Introducing the Burgers form b = & a (x)-$ a = b^xjdx 11 , together with Eq. [21], 
this may be rewritten as 

b = -Q + tt L . (35) 
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This is a very important equation. It relates the bond angle, computed 
from the actual connection Q as indicated in Eq. |25], to the distribution of 
dislocations and to the Gaussian curvature of Eq. |33|. We can make this 
statement more apparent by applying the exterior derivative to Eq. [35], and 
using Eqs. |29| and [33|: 

K{x) = s(x) + -^^(x) (36) 



o 1 N 01 

-K(*) = ^E^(x-x,) + -^e^yax). (37) 
71 yj\9\i=i \J\9\ 

Note that the coordination number of an atom n,, the number of nearest- 
neighbors, is given by 

f 1 N 3 f 3 f 
Q-n i = q i= / du^/g—i= ^ g^(x - x,-) = - / du^fgK / b^dx^ , 

JY1 . q j 7T jE 7T J OS 

(38) 

where Sj is a small region around the atom. This result contains an additional 



surface term when compared with the existing literature Wu, 11]. This surface 



term will only vanish if defects successfully fully screen out the Gaussian 
curvature. As we will see in Eq. 59, this is the condition that minimizes the 
energy of the model. So, at very low temperatures, the surface term will be 
small. 

Comparing Eq. [37] with [jllj, there is a somewhat disturbing factor ^ in the 
last term. In a fine grained description of the model, one regards a dislocation 
as a composite object, a tightly bound pair of opposite disclinations (see 
Fig. [5p. Defining the dipole moment p as 

= l^p , (39) 

7T 

where b v is the microscopic Burgers vector, we can Taylor expand the density 
of two free disclinations of opposite charge as 

I {5(x + p) - 5(x)} = e^fi^x) + 0(a 2 ) , (40) 

where we neglect higher order terms in the lattice spacing a. The right hand 
side of this last equation corresponds exactly to an isolated dislocation. The 
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factor - is the slight displacement of the center of the disclinations necessary 
to render the right dislocation density Eq. [4(| As a byproduct, we derive the 
result that the dipole moment p is perpendicular to the microscopic Burgers 
vector, a result which is also apparent from Fig. |3|. 




Figure 3: A dislocation as a tightly bound disclination pair. The dot- 
ted line is the Burgers circuit. It is apparent that the dipole moment p is 
perpendicular to the Burgers vector b. 

For our purposes, however, it will be more convenient to use Eq. |35] in 
combination with Eq. |31], which yields 

b = -de- n sing + n L . (4i) 

This completes our geometrical development of the structure of bond- 
orientational order. 

4 The properties of the hexatic phase 

The previous section was somewhat formal. We derived some general results 
valid in the presence of bond-orientational order. Further progress requires 
making additional physical assumptions about the interactions of dislocations 
and disclinations with the underlying geometry and topology. 

For the simplest case of a flat monolayer a free energy was proposed 
in m |2| (see for a review). We now examine this proposal within the 
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framework of the geometric formalism developed here. We use a completely 
covariant language with the aid of the mathematical concepts sketched in 
Sec. [| Although this is not necessary for the case of flat space, it is directly 
applicable to general geometries. 

4.1 Debye-Huckel approximation on a flat monolayer 

On a flat monolayer, the metric of Eq. p] is the standard Euclidean metric 

ds 2 = dx 2 + dy 2 . (42) 

In this case i? 1 = dx and d 2 = dy. 

Let us first treat the simpler case of vanishing free disclination density. 
The bond angle as a function of the Burgers vector density is found by 
inverting Eq. |4l]: 

e^^^-Lj^J^l, (43) 

where \ = 1/(2tt) log() x — x'\). This result is in agreement with the results 
in |L], 0, |3|. We emphasize here that this this result follows from the Cartan 
structure equations alone, with no further physical assumptions. It provides a 
non-trivial consistency check for our geometrical interpretation of the hexatic 
phase. 

Additional constraints follow from Eq. With no free disclinations it 
implies 

db = t^d^K = . (44) 

Geometry thus constrains the Burgers vector distribution b(x) to be irrota- 
tional. 

The free energy of the crystal, within a Debye-Huckel approximation, was 
derived in Jj], @: 

n/T=^(db\-^db) + ^(b\b) , (45) 

where the inner product between forms was defined in Eq. |6] and the Laplacian 
in Eq. Kb is the 2-dimensional Young modulus, a is the original lattice 
spacing and E c the core energy of a single dislocation. The RG recursion 
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relations of the 2d KTNHY melting transition []J imply that in the hexatic 
phase 

\im p ^ K B (p) = . (46) 

If now we consider the bond angle as the fundamental variable, Eq. [H| and 
Eq. ^4| imply a simple XY model for the energy of the bond angle variable 

n/T = ^J d^e , (47) 

with hexatic stiffness K A = ||, §. 

The irrotationality of b is a general constraint to be satisfied in any 
hexatic regime with vanishing free disclination density. It has additional 
consequences. We already noted that a dislocation may be regarded as a 
tightly bound disclination pair. The dipole vector p, defined in Eq. [I!], is 
perpendicular to the Burgers vector b. If we consider a closed loop C, we 
can compute the circulation for p 

J pds = J e^budxp = J dxdyAO , (48) 

where S is the two dimensional region enclosed by the closed curve C. The 
resultant integral is generally not zero. This means that dipoles form grain 
boundaries of consecutive fives and sevens. If there are no free disclinations 
these strings cannot self- intersect. 

In the presence of free disclinations, Eq. gets modified to 

db = sdx A dy e^d^b v {x) = s(x) . (49) 

The free energy for the bond angle variable now reads 

h/t = ^(s\^s) + ^(de + n smg \de + n sing ) . (50) 

The first term leads to a strong binding of free disclinations, since it corre- 
sponds to a |x| 2 ln |x| disclination-disclination interaction . Recalling Eq. |46| , 
however, we see that this term is absent. The long-wavelength description of 
the hexatic phase is given by the standard XY model ffl]. 



The hexatic phase can be characterized by the correlator [12 



C(v) = (a(v)a(0)) , (51) 
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with 

<7(r) = Ind /^( r " = ^( x ) , (52) 
where Ind/ is the index of the order parameter. The last identity follows 



from the definition of the index of a vector field. The order parameter Eq. 51 



is thus the correlation function of the free disclination density, the hexatic 
curvature. To lowest order in the fugacity, we have in the hexatic phase 

< S ( X M°)> ~ > (53) 



where y is the fugacity of the vortices. In the isotropic liquid phase [12 



(s(x)s(O)) ~ e~~, (54) 

a 

where £ ~ e T ~ T c is the correlation length of the model. 



4.2 Hexatic order in an arbitrary geometry 

In the previous section we recast standard results for a flat monolayer in the 
geometric framework developed in this paper. These results may now be 
extended directly to the case of an arbitrary geometry. The metric 

ds 2 = g^dx^dx" , (55) 

takes a general form with metric coefficients given by Eq. [TJ. 

The first step is to find the free energy. This is easily obtained from 
Eq. |45| using the covariant definition of the operators involved, as described 
in sect. |2] We start by defining the form 

p(x) = {K(x)-s(x)}fi M , (56) 

expressing the difference between Gaussian curvature and free disclination 
density. In Eq. |45| we deliberately wrote the free energy of a simple flat 
monolayer in terms of the torsion degrees of freedom. It is now extremely 
simple to generalize it to the case of a general geometry by making use the 
general relation Eq. |5l] that holds for an arbitrary geometry to obtain 

H/T = *±(p\^Lp) + ^(d9 + Q sing - Q L \d9 + Q sing - Q L ) , (57) 



15 



where Vl L is the connection form of the Levi-Civita connection. This corre- 
sponds to the free energy of a 2D crystal on an arbitrary geometry. This 
result is identical to that obtained by integrating out the in-plane phonons 
in the absence of free disclinations, and generalizing to include these ad- 
ditional degrees of freedom. The free energy Eq. |57| corresponds interacting 
dislocations and disclinations in the crystalline phase. At zero temperature, 
one can substitute the bond orientational order parameter 9 by its minimum 
value 9 = ^V^f^, and we get 

H/T= x (pl i p) + : r (pl i p) - (58) 

Free positive disclinations are attracted to positive curvature regions and 
negative disclinations to negative curvature regions. The ground state is 
defined by the equation 

p = s (x) = K{x) . (59) 

Defects will arrange themselves to optimally screen out the Gaussian curva- 
ture. 

As for the case of a flat monolayer (Eq. [|6|), the hexatic phase requires 
that in the infrared limit 

K B = . (60) 

The free energy then becomes 



T~L/T = — ^-(d9 + VL S i ng — Vt L \d9 + VL sing — Vt L 



K_a, 
2 

^ J dxy/ggV&J + n sing - fij)(0„0 + Vising - V L U ) , (61) 



which is the hexatic free energy first considered in ||. 

It is beyond the scope of this paper to investigate the precise mechanism 
by which Eq. ^0| is implemented in an arbitrary geometry, but we show instead 
that it is a necessary condition for the hexatic phase to be realized. 

In analogy with the flat monolayer, one may compute the order parameter 



Eq. 51 to lowest order in the fugacity, 



(s(x)s(O)) = -2(|)%V^ Exi ^ (i) ^ie^SL 1 /^'f^^(x) ; 

(62) 
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which reduces to Eq. ^ in the case of a flat monolayer. 

It is worthwhile to recall at this point that if fluctuations in the geometry 
were allowed, we should include an additional bending rigidity or extrinsic 
curvature term 

^ = \ J dx^gH* , (63) 

where H is the Mean curvature of the surface. In this case we average over 
all possible fluctuations in geometry. It seems that Dislocations have a finite 
energy || and therefore they proliferate at any temperature, see also 0|. 



5 Conclusions and Outlook 

A full understanding of defects in curved geometries is a challenging subject 
with many novel features not encountered in the plane. In this paper, we 
presented a geometrical interpretation of bond-orientational order leading to 
a variety of relations between the underlying geometry, the topology and the 
nature of the defects. The formalism developed in this paper has already been 
applied very successfully to the case of defect arrays in spherical crystals fll? 



Further applications are currently being explored. The most challenging open 
problem may be the generalization of the KTNHY renormalization group 
flows to fixed or dynamical curved geometries. The spherical crystal itself 
gives rise to a very rich structure of defects in the ground state |T7| and other 
geometries are very likely to lead to even more novel results. 
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A The hexatic cone 



A sample in the hexatic phase on a template having the geometry of a cone 
is the simplest case where the present formalism may be applied. It provides 
an interesting benchmark as the cone has been the subject of much interest 



1% |15| , |16| . For some range of parameters, an isolated disclination forces 
the flat monolayer to buckle; a positive disclination is optimally screened by 
the Gaussian curvature located at the cone. It is more convenient to work in 
polar coordinates (r,ip), the metric being 

ds 2 = (1 + m 2 )dr 2 + r 2 dip 2 , (64) 

where m is related to the aperture angle of the cone ( by tan£ = — . 

One readily finds i? 1 = a/1 + m?dr, § 2 = rdijj. The Levi-Civita connection 
is Q L = — v / 1 ^ m 2 dip. The actual hexatic connection is flat, corresponding to 
a plane 

Q= -dif> + d<l>^>e(r,ip) = -ip + <f){r,ip) , (65) 

where is regular 0(r, ip + 2tt) = 0(r, ip) and corresponds to the different 
Burgers vectors distributions one may have in a cone. The lowest energy 
solution is = 0, and this is the only case we consider. The torsion vector 
distribution may be computed easily, 

b(x) = (l- 7 =L=);e v ,, (66) 

Torsion vectors are tangent to circles centered at the tip of the cone, as 
depicted in Fig. f|. It is also illuminating to plot the order parameter, as 
illustrated in Fig. ^. The parallel transported bond angle is completely in- 
sensitive to the Gaussian curvature located at the tip of the cone. As a further 
cross check, we can compute the hexatic energy for this configuration, using 
the free energy Eq. ^l[ Assuming a lattice spacing a as an ultraviolet cut-off 
and the radius R of the cone as an infrared cut-off, we get 

Ka f R dr , 9x i f 2n , . 1 n9 

E hex = -± -(l + m 2 W ^(-l + ^== 2 
2 J a r Jo v 1 + m 2 

= ^(l--=^=) 2 (l + m 2 )^ ln ^. ( 67 ) 



This result has already been obtained by other authors [15], [L6| by solving 
the equations of motion for the order parameter and computing the Green's 
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function. Instead we have obtained it by constructing the order parameter 
by parallel transport. As a byproduct, we are able to compute the torsion 
vector distribution in the cone. 




Figure 4: Distribution of torsion vectors in the cone. Plot of the order 
parameter in polar coordinates (r, ip) . The order parameter is the same as in 
a plane despite the Gaussian curvature located at the tip T of the cone. 

If there is an isolated disclination located at the tip of the cone, the 
hexatic connection is no longer flat, as there is hexatic curvature located at 
the tip of the cone. Parallel transport is now ambiguous if we encircle the 
tip of the cone. We get, 

Q = (§ - 1)# -> 9(r, = (f - lty . (68) 
o o 

Following the same steps as before we obtain a free energy of the form given 
in Eq. 
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